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A novel technique to measure heat transfer and liquid film thickness distributions over relatively large 
areas for two-phase flow and heat transfer phenomena using infrared (IR) thermometry is described. 
IR thermometry is an established technology that can be used to measure temperatures when optical 
access to the surface is available in the wavelengths of interest. In this work, a midwave IR camera 
(3.6-5. 1 pm) is used to determine the temperature distribution within a multilayer consisting of a silicon 
substrate coated with a thin insulator. Since silicon is largely transparent to IR radiation, the temperature 
of the inner and outer walls of the multilayer can be measured by coating selected areas with a thin, IR 
opaque film. If the fluid used is also partially transparent to IR, the flow can be visualized and the liquid 
film thickness can be measured. The theoretical basis for the technique is given along with a description 
of the test apparatus and data reduction procedure. The technique is demonstrated by determining the 
heat transfer coefficient distributions produced by droplet evaporation and flow boiling heat transfer. 

© 201 1 Elsevier Ltd. All rights reserved. 


1. Introduction 

Definitive understanding of phase change heat transfer mecha- 
nisms remains elusive due its sensitivity to many variables, but also 
due to a lack of reliable local information that can enable models to 
be tested. Although point measurements or area averaged measure- 
ments of variables such as local film thickness and heat transfer have 
been made, techniques whereby these quantities can be measured 
over large areas are generally lacking. Point or average measure- 
ments may be appropriate at very low heat fluxes where insignifi- 
cant evaporation occurs and the flow regime does not change, but 
such measurements are insufficient for model validation. 

Techniques that have been used in the past to measure local 
film thickness have utilized capacitance sensors, conductance 
probes, confocal microscopes, and other techniques such as reflec- 
tance (e.g., Coney, 1973, Klausner et al., 1990, Han and Shikazono, 
2009, and Shedd and Newell, 1998). The heat transfer has usually 
been measured using thermocouples welded to the tube walls or 
by resistively heating the walls and measuring the average wall 
temperature. The authors are not aware of any techniques whereby 
liquid film and heat transfer distributions in flow boiling are mea- 
sured over relatively large areas with high resolution. If the heat 
transfer distribution along the walls of a tube could be measured, 
it could be used to verify models of the heat transfer variations pre- 
dicted by models of slug flow, wavy annular flow, annular flow 
dryout, etc. For example, a three-zone model of elongated bubble 
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evaporation in microchannels was proposed by Thome et al. 
(2004) and Dupont et al. (2004). The model assumes that heat is 
transferred to a liquid slug, an elongated bubble, and a vapor slug 
(Fig. 1 ). Heat transfer to the elongated bubble is thought to occur 
through the thin liquid film between the vapor and the wall. The 
bubble frequency, liquid film thickness, and heat transfer to the li- 
quid and vapor slugs are obtained from correlations. The globally 
averaged heat flux was measured and used for model verification. 
Much stronger verification of the model can be made if measure- 
ments of the local heat flux were available to verify the heat trans- 
fer rates before, during, and after passage of the bubble. 

In flow regime based models (e.g., Kattan et al., 1998a,b,c), the 
local flow boiling heat transfer coefficients are predicted based 
on the tube perimeter fraction wetted by liquid. For example, for 
a horizontal stratified flow in a circular tube, the local heat transfer 
coefficient is given by 

i E id dry ^l vapor + r x (2 7T — 9dry)h we t , 


where h wet is assumed to be composed of a nucleate boiling compo- 
nent and a convective component in a power law form: 



For annular and intermittent flow, h tp — h we t • The local void frac- 
tion needed to evaluate liquid and vapor velocities is calculated 
using correlations (e.g., Rouhani and Axelsson, 1970). Again, verifi- 
cation of these models can benefit greatly from measurements of 
the local heat transfer distribution in the various regimes. 
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Fig. 1 . Three-zone heat transfer model for evaporation in a microchannel from [1 ]. L p is the total length of the triplet. L h L dry and are the length of the liquid slug, the dry 
zone and the liquid film trapped by the bubble, respectively. The internal radius and diameter of the tube are R and d while < 5 0 and S min are the thicknesses of the liquid film 
trapped between the elongated bubble and the channel wall at its formation and at dry out. 


The objectives of this paper are to describe and demonstrate an 
IR thermometry based technique whereby the local heat transfer 
produced by two-phase flows can be measured with high spatial 
and temporal resolution. The application of the technique to the 
measurement of liquid film thickness is also described. IR ther- 
mometry is an established technology that can be used to measure 
temperatures when optical access to the surface is available in the 
wavelengths of interest, and has been used to measure heat trans- 
fer distributions during pool boiling heat transfer. Golobic et al. 
(2009) and Stephan and co-workers (e.g., Schweizer and Stephan, 
2009) used an IR camera to measure the heat transfer distribution 
under single nucleating bubbles as they grew on a thin metal foils. 
Stephan has recently begun using a thicker CaF substrate in place 
of the thin film in order to increase the heat capacity of the sub- 
strate (Fischer et al., 2011) so it is more representative of real sur- 
faces. Gerardi et al. (2010) used a high speed IR camera in 
conjunction with a video camera to measure bubble behavior on 
an ITO heated sapphire substrate. The IR camera measured the 
temperature distribution at the ITO surface, while a video camera 
was used to visualize the fluid behavior. Krebs et al. (2010), Shen 
et al. (2010) and Mani et al. (2012) used IR thermography to study 
flow boiling in microchannels, droplet evaporation, and jet 
impingement, respectively. In these studies, an IR camera was used 
to view through a silicon substrate to visualize in great detail the 
temperature distribution at the silicon/water interface. Sefiane 
et al. (2008) used IR thermography to visualize the spontaneously 
occurring hydrothermal waves within evaporating methanol, etha- 
nol, and FC-72 droplets. 

In the present work, a midwave IR camera (3.6-5. 1 pm) is used 
to measure the temperature variations within a multilayer consist- 
ing of a silicon substrate coated with a thin thermal insulator that 
is partially transparent to IR. The insulator amplifies the tempera- 
ture variations and provides a strong signal for the IR camera. Since 
silicon is largely transparent to IR radiation, the temperature of the 
inner and outer walls of the multilayer can be measured by coating 
selected areas with a thin IR opaque film. The fluid used (FC-72) is 
also partially transparent to IR over a broad range of wavelengths, 
allowing the flow to be visualized and the film thickness to be 
measured. The theoretical basis for the technique, a description 
of the test apparatus and data reduction procedure, and experi- 
mental validation are presented in the sections below. 

2. Theoretical background 

Consider the multilayer wall consisting of a silicon substrate 

onto which polyimide tape (polyimide layer + acrylic adhesive) is 


attached as shown schematically on Fig. 2a. An opaque black paint 
that is much thinner than the other layers is applied to the top of 
the polyimide tape. The black surface is exposed to a two-phase 
flow. The polyimide tape is necessary to measure heat transfer 
coefficient distributions of the expected magnitude since the high 
thermal conductivity of the silicon substrate would simply smear 
out any temperature variations through substrate conduction, 
reducing both the magnitude of the temperature differences as 
well as the spatial resolution. 

To obtain the heat transfer coefficient at the fluid-wall interface, 
the temperature gradient within the polyimide tape is required. If 
the time-varying temperatures of the black surface T s i(t) and T s2 (t) 
are known, the instantaneous temperature profile within the mul- 
tilayer can be obtained through an unsteady heat conduction sim- 
ulation. If the polyimide tape is thin compared with the spatial 
resolution of the camera and the temperature gradient is much lar- 
ger in the x-direction than in the y- and z-directions, a 1-D heat 
conduction assumption can be used. The 3-D heat conduction 
equation must be used otherwise. Assuming 1-D heat conduction, 
the governing equations within the layers are given by 

Psi^p si T^T — ksiV T + Qsi (3a) 

PaCpA^ = '< a V 2 * T (3b) 

PpCp,p'-^ = kpS7 2 T (3c) 

and the system is subject to the boundary conditions T=T s i(y,z,t) at 
x = 0,T= T s2 (y,Z, t) at x = L si + L A + L P . 

Consider now the calculation of the black surface temperature 
T sl (t). This temperature is not directly available since the energy 
measured by the IR camera consists of emission from the black sur- 
face, emission from each of the layers (which depends on the tem- 
perature profile within them), and reflection from the 
surroundings. Since the optical properties of the polyimide and 
adhesive are similar, they are treated as a single layer in the radi- 
ation calculation as indicated in Fig. 2b. The total energy measured 
by the camera (E c ) is the sum of the energies emitted by each layer 
within the spectral bandwidth of the IR camera (2 i-2 2 ): 

Ec = Poo-c^oo ~"F Osi-c^si 6j-cEj + T s-cEs (4) 

where E ^ _ k2 <jT^ is the blackbody radiation due to the sur- 

roundings, E si = / 0 4 ccsiFx,-x 2 <J[Tsi(,x)f exp( y. si x)dx is the energy 
emitted by the silicon that reaches the Si-oo interface, 
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Fig. 2. Model descriptions, (a) System schematic for the conduction problem, (b) System schematic for the radiation problem. 


E T = Jq T oc t F^-a 2 (t[Tt(x)} 4 exp(-ot T x)dx is the energy emitted by the 
tape that reaches the T-Si interface, E s = F^-^crT^ is the blackbody 
radiation of the black surface. 

F^~x 2 is total emitted energy from a blackbody contained with 
the wavelength interval X\-X 2 and can be obtained from tables or 
by integrating the Plank distribution. The coefficients p^, s si _ c , 
s T -c, and t s _ c account for the attenuation and reflection within the 
multilayer wall, and are given by 


P 00 — C = 

0 Psi-oo ) Papp.Si-T^ Si 
~ Psi-oo + 1 _ n n r 2 

1 Psi-oo Papp,Si-T L Si 

(5) 

^Si-c = 

(1 — Psi-oo)0- Papp.Si-T^'Si) 

(6) 

1 — Psi-oo PappSi-T^Si 

£r-c = 

(1 — Psi- oo)(l — Psi-T )(1 + Pr-s^rj'Fsi 

(7) 

(l - Psi-ooPapp,Si-T T Si) 0 “ Psi-rPr-stf) 


Ts-c — 


0 ~ Psi-JO - Psi-t)0 - Pj^siTr 
(1 — Psi-ooPapp^i-T^Sij 0 — Psi-rPr-stf) 


Papp y Si—T — Psi T + 


(1 ~ Psi-r) Pt s T T 
1 ~ P Si— T P t—s^T 


( 8 ) 


(9) 


The derivation of the above equations is given in the Appendix. 
The optical properties of the various layers along with the temper- 
ature of the surroundings (T^) are assumed to be known. The tem- 
perature distribution within the silicon [T Sl (x)] and the tape [Tt{x)] 
and the temperature of the black surface (T s> i) are not known ini- 
tially, but can be obtained by solving the coupled conduction and 
radiation problems according to the following algorithm: 


1. Assume an arbitrary temperature profile within the multi- 
layer at t= 0. 

2. Compute E S i and E T from the temperature distribution and 
determine an updated E s and surface temperature T s , i from 
Eq. (4). 

3. Solve the conduction equation using the updated T St i to 
obtain a new temperature profile at t = At. 

4. Repeat steps 2 and 3 for each successive time step. 

The effect of the assumed initial temperature profile within the 
multilayer will decay after which the true temperature profile will 
be known. 1 The heat flux from the wall to the fluid can be obtained 
from the derivative of the instantaneous temperature profile within 
the polyimide at the black surface according to q" = k P |£| . 

3. Experimental results 

3.1. Radiant property determination 

The above model requires knowledge of the reflectivity and 
absorptivity of the materials used. These were determined experi- 
mentally using the two configurations shown in Fig. 3. The black- 
body consisted of a large cylindrical cavity at a controlled 
temperature behind a small orifice as shown in Fig. 4. The emissiv- 
ity of this blackbody is estimated to be >0.999 using the relation 
provided in Quinn (1967). 

The apparent reflectivity was determined using the experimen- 
tal setup shown in Fig. 3a. The blackbody was positioned such that 
emission from it reflected from either the front or back surface of 
the layer to be tested. The total energy reaching the IR camera con- 
sisted of the blackbody energy reflected from the layer, the self- 


1 For the condition used for the experimental verification described below, the 
initial transient was found to have completely decayed after 0.1 s. 
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emission of the layer, and the transmission of energy from the sur- 
roundings through the layer: 

Eci — Papp,m — co (^i ~h H - £app.m — oo (j* h~?-2 

+ Tapp,m- 00 (F^_ X2 aT 4 x ) ( 10 ) 

The measurement was then repeated without the blackbody: 

£c 2 = Papp. f m- 00^1-^2^00) + e app,m-oo^i_^ 2 CtT^ 

+ ^app : jn- 00 -;. 2 ) 0 ^ ) 


Subtraction of Eq. ( 1 0) from Eq. ( 1 1 ) allowed the apparent reflectiv- 
ity to be determined since the temperature T b and T ^ are known. 

The apparent transmissivity was determined using the setup 
shown in Fig. 3b. The total energy reaching the camera with and 
without the blackbody are given by 



Subtracting one from the other allowed the apparent transmissivity 
to be determined. The apparent reflectivity and transmissivity for a 
single layer m sandwiched by two identical layers n are given by 


, + 


0 Pm- 00) Prr 
WLoo? 


(14) 


T 


app..m - 00 


(1 ~ Pm- 00) 

1 - o 2 T 2 

1 rm-oo L m 


(15) 


(see Appendix for a derivation), from which p m -oo and r m can be ob- 
tained. The absorption coefficient can then be obtained from the 
Beer-Lambert law, T m , and the layer thickness. 

The above technique was verified by measuring the reflectivity 
of a double side-polished silicon wafer (n = 3.43 at k = 3.7 pm) in 
air (n = 1.00). An analytical expression for the reflectivity at the 
silicon-air interface is given by 


Psi-air 


f nsi - n Qi r \ 2 

\n S i + n air ) 


0.301 


(16) 


This was in good agreement with the experimentally measured va- 
lue Of psi-air =0.34. 

The reflectivity of the tape was calculated using the indices of 
refraction of the polyimide {n po iyimide = 1-7) and air to be 



Fig. 4. Blackbody construction. 


P polyimide- air 


f ft polyimide H-air^ 

yl polyimide ~\~ Hair, 


0.07 


(17) 


If the index of refraction of the acrylic adhesive was assumed to 
be that of acrylic (n a dhesive = 1.5), the reflection would be 
P adhesive-air = 0.04. The reflection between the polyimide and adhe- 
sive is very small p po iyimide-adhesive = 0.004. Since reflection within 
the tape is small and since the reflectivities of the polyimide and 
the acrylic adhesive are similar, the tape was considered to behave 
as a single optical layer with polyimide properties. The absorption 
coefficient of the tape was measured experimentally. A summary 
of the optical properties of the silicon and polyimide tape are given 
on Table 1. The absorption coefficient of the polyimide tape is seen 
to be much larger than that of the silicon. 


3.2. Uncertainty analysis 

The experimental uncertainty in the calculated heat flux was 
determined by perturbing each experimental parameter a small 
amount one at a time and running the algorithm to find the sensi- 
tivity to the variables at varying heat fluxes. Many of the assumed 
physical constants (e.g. thermal conductivity and thickness of sili- 
con) introduced errors that were several orders of magnitude smal- 
ler than the important sources of error and were therefore omitted. 
The parameters that resulted in the highest error along with their 
values and assumed are tabulated in Table 2. The maximum total 
error, calculated as the RMS of each individual error, increases with 
heat flux and is summarized in Table 3 at various heat fluxes. 


3.3. Validation and demonstration 
3.3.1. Air cooling 

To validate the experimental technique and algorithm, the 
apparatus shown in Fig. 5 was used. A Si wafer was heated by a 
thin film heater. Polyimide tape (30 pm total thickness) with a thin 
black coating was deposited onto the top of the Si wafer to 


(a) 


Film, T m 


(b) 


Film, 


Blackbody 


1 E„ 


"~1r h +— £ - 


IR 


e „ r 

L 


Blackbody 


Fig. 3. Schematic of experimental setup used to measure the reflectivity and transmissivty of a single layer, (a) Apparent reflectivity measurement, (b) Apparent 
transmissivity measurement. 
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Table 1 

Optical properties of the silicon and polyimide tape. 


Material 

Absorption coefficient (m 1 ) 

Silicon 

52.6 

Polyimide tape 

7110 

Interface 

Reflectivity 

Si-air 

0.34 

Polyimide-air 

0.07 

Polyimide-Si 

0.14 


measure the top surface temperature I s i(t). The black coating 
(6 pm thick) was produced by silk screening onto the polyimide 
tape Nazdar GV111 paint consisting of about 20% carbon black in 
a vinyl chloride/vinyl acetate copolymer after curing. The high car- 
bon black content ensured the coating was opaque, and its high 
thermal conductivity relative to the polyimide and thinness al- 
lowed its temperature to be assumed to be uniform. 

The emissivity of the black surface was experimentally deter- 
mined. The polyimide tape with black coating was attached to a 
heated aluminum plate containing a thermocouple. The plate 
was heated to temperatures between 45-95 °C, and the energy 
emitted from the surface along with the energy reflected from 
the surroundings were measured using the IR camera. The total en- 
ergy measured by the IR camera is given by 

E c = sF h aT 4 s + (1 - s)F h gT 4 x (1 8 ) 

from which the surface emissivity can be determined. The emissiv- 
ity of the black coating was measured to be 0.90. 

A small piece of a polyimide tape with black coating was at- 
tached to the bottom of the wafer to measure the bottom surface 
temperature T s2 (t). Gold coated mirrors were used to simulta- 
neously measure the energy emitted from both the top [E Top (t)] 
and bottom [E BottC {t\E Bottb {t)] surfaces. Since the black surfaces 
were not transparent to IR, the black surface temperatures on the 
top [T sl (t)] and bottom [I s2 (t)] surfaces could be directly obtained 
using the energies [E Top {t\E Bot)b {t)] using an IR camera (Electro- 
physics Silver 660 M). The black surface temperature on the top 
of the wafer T sl (t) could also be calculated based on the algorithm 
using E Bott c{t) and T s2 (t). To validate this calculation, the calculated 
top surface temperature was compared with the temperature ob- 
tained from E Top (t) as the top of the wafer was cooled using an 
air jet. Results are shown in Fig. 6. The air jet was turned on at 
t = 5.7 s, resulting in a drop in the surface temperature. When the 
air jet was turned off at t = 10.9 s, the temperature recovered. Good 
agreement within 0.5 °C was seen between the measured and com- 
puted values, demonstrating that the algorithm can be used to 
accurately predict the true surface temperature. 

3.3.2. Droplet evaporation 

To further demonstrate the power of the technique, the appara- 
tus shown in Fig. 5 was used to measure the wall heat flux during 


Table 2 

Parameters with largest effect on uncertainty analysis. 


Experimental constant 

Value 

Uncertainty 

Emissivity of the black coating 

0.90 

0.01 

Absorptivity of the polyimide tape 

7110 (m- 1 ) 

500 (m- 1 ) 

Polyimide thickness 

15 pm 

2 pm 

Adhesive thickness 

15 pm 

2 pm 

Reflectivity of silicon-air interface 

0.34 

0.017 

Reflectivity of silicon-polyimide interface 

0.14 

0.007 

Thermal conductivity of polyimide 

0.12 W/m-K 

0.01 W/m-K 

Thermal conductivity of adhesive 

0.20 W/m-K 

0.01 W/m-K 


Table 3 

Summary of uncertainty analysis. 

Heat flux (W/cm 2 ) 2.0 4.0 6.0 8.0 10.0 

Max Error (W/cm 2 ) 0.35 0.47 0.62 0.80 0.98 



Fig. 5. Experimental setup for validation (not to scale). 



Fig. 6. Comparison of the top temperature, Tsi(t), determined using the algorithm 
and the energy measured through the bottom £ BotiC (t ) (computed temperature) with 
that measured directly using E top (t ) (actual temperature). 

evaporation of a single droplet of PF-5060. The substrate was 
heated to 83 °C by attaching a thin film heater to the silicon wafer, 
after which a syringe was used to deposit a droplet onto the top 
surface. The IR camera was used to measure the heat transfer to 
the droplet from below the wafer. Fig. 7 shows temperature and 
heat flux contours as the droplet evaporates and illustrates how 
the local temperature and heat transfer can be measured with high 
spatial and temporal resolution. 2 

The droplet reached its first maximum diameter approximately 
21 ms after impact, then contracted until 42 ms. During this pro- 
cess, high heat transfer between the droplet and the heated surface 
was observed inside the droplet with a maximum heat transfer of 


2 A movie of the droplet evaporation process is available online at (mmcl.mpg). 



T.H. Kim et al./ International Journal of Multiphase Flow 40 (2012) 56-67 


61 



Temperature (°C) 








t = 273 ms 


Temperature (°C) 



Heat Flux (W/cm 2 ) 



t = 308 ms 



t = 317 ms 



t = 329 ms 



t = 336 ms 



t = 368 ms 



t = 399 ms 



t = 459 ms 


Fig. 7. Temperature and heat flux contours for droplet evaporation. 


about 23 W/cm 2 . After 42 ms, the droplet spread again and the 
diameter reached a second maximum at 242 ms. During this time, 
nucleate boiling within the droplet resulted in a higher heat trans- 
fer (about 10 W/cm 2 ) at the center of the droplet than at the edges. 
After 242 ms and until the droplet evaporates, the droplet size de- 
creased due to evaporation at the rim of the droplet. Nucleate boil- 
ing ceased as the droplet thinned and ruptured. The highest heat 
transfer occurred at the receding three-phase contact line as 
expected. 

The above results are for PF-5060 evaporating on the black 
painted surface and may not be representative for evaporation 
on silicon. Due to the highly wetting nature of PF-5060, however, 
the effect of contact angle is not expected to be important. Surface 
effects may be much larger when water is used, however. 


3.3.3. Gravity effects on flow boiling 

The technique was also used to measure flow boiling heat trans- 
fer during the high gravity and low gravity conditions produced by 
an aircraft. A polished tube 6 mm ID x 8 mm OD x 120 mm long 
was manufactured from a block of single crystal n-type silicon with 
>1 Q-cm resistivity. Two current taps and two voltage taps consist- 
ing of a chromium adhesion layer and 2000 A of aluminum were 
deposited at the ends of the tube after they were doped to mini- 
mize contact resistance (<1 Q), allowing the power dissipated by 
the tube to be determined using a four-wire measurement. The 
tube was heated using a power supply that could impose up to 
300 V across the tube, producing a nominal heat flux of up to 
4 W/cm 2 . Polyimide tape was attached to the inside of the tube, 
but only half of the tape had a black coating (Fig. 8). The mirror 
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Fig. 8. Mirror setup to allow simultaneous flow visualization and measurement of 
multilayer temperatures. The green lines indicate the location of the black coating 
on the inner and outer surfaces of the silicon tube. (For interpretation of the 
references to color in this figure legend, the reader is referred to the web version of 
this article.) 

arrangement allowed simultaneous visualization of the flow along 
with measurement of the temperature distribution on the inside 
and outside of the tube. 

A flow loop (Fig. 9) was constructed so FC-72 could be pumped 
through the tube at known mass fluxes up to 200 kg/m 2 -s and inlet 
temperatures from room temperature to near saturated 
( T S at = 57 °C for FC-72 at 1 atm). The primary loop included the sil- 
icon test section where wall temperature measurements were 
made using the IR camera. The two-phase mixture exiting the test 
section was condensed in a counterflow heat exchanger then 
pumped through a flowmeter and preheater before entering the 
test section again. An accumulator (expansion tank) was used to 
keep the pressure within the loop at nominally 1 atm. The temper- 
ature of the FC-72 exiting the heat exchanger was controlled by 
occasionally pumping cold water through the heat exchanger. 

The variable gravity environment was produced by a Zero G 
Corporation 727 aircraft flying parabolic arcs as part of NASA’s 
FAST program in July, 2011. Approximately 20 s of 1.8 g pullup 
was followed by 20 s of low-g (±0.01 g) and 20 s of 1.8 g pullout. 
Up to 40 parabolas were flown per flight and four flights were 
made. The test section was oriented such that the liquid entered 
the tube in a vertical upward flow configuration— because the flow 
produced is nominally axisymmetric, heat transfer data was re- 
duced only along the axis of the tube. This configuration also 
avoids refraction effects that can occur due to the curvature of 
the tube when viewing the tube off-axis. 

The heat transfer coefficient distribution is shown on Fig. 10a 
as the fluid flows through the tube in a low-g environment at a 
relatively low mass flux of 44 kg/m 2 -s. 3 Large variations are ob- 
served where liquid rewets the surface with peaks up to 
12,000 W/m 2 -I<. Very low heat transfer is also observed on the 
tube, indicating dryout of the surface. Time averaged heat transfer 
coefficient distributions along the tube in high-g and low-g are 
shown on Fig. 10b. The effect of gravity is minimal in the first half 
of the tube where the flow is dominated by bubbles dispersed in a 
liquid. The heat transfer deteriorates towards the tube exit, how- 
ever, as the flow transitions to a wavy annular type flow with peri- 
odic dryout. Additional data at various flow rates and inlet 
enthalpies will be presented in future work. 


3.3.4. Liquid film thickness 

We close this paper with a discussion on how the technique can 
be used to measure liquid film thickness in a tube. Consider the 
case where a liquid film flows along the inside walls of a round 


3 A movie of flow boiling in low-g is available online at (mmc2.mp4). 



Test Section 


Gear Pump 


Water/Ice Box 


Preheater 


Heat Exchanger 


Water Pump 


Fig. 9. Solid model of flow loop. 


tube as can occur during annular flow (Fig. 11). The tube and liquid 
are in thermal equilibrium at temperature T s , and placed in a large 
enclosure at temperature T^. The IR camera looks through the tube 
at a black surface at T hot - The energy measured by the IR camera E c 
consists of energy from the surroundings that are reflected from 
the front, the emission from the tube and liquid, and the energy 
emitted from T hot that is transmitted through the tube and liquid. 
It can be shown that 

Ec = Poof— C Eb } oo + [(1 — Tm)£m/-c H"~ (1 — T l)fyf-c H"~ (1 — T/)S/r-c 

+ (1 — 't m )£mr-c]Eb > S T" / Lhot/-cEb..hot ( 19 ) 

where 

E b ,oo = <ylf, E bi s = aT 4 s , E b>hot = oT\ 0t 
t, = exp(-a,L,), T m = exp(-a m L m ) 
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Time = 35.475 sec 




Fig. 10. Flow boiling heat transfer during vertical upward flow at a mass flux of 44 kg/m 2 -s and fluid inlet temperature of 55 °C. The heat flux applied to the tube is nominally 
4 W/cm 2 . (a) Heat transfer coefficient distribution at a representative time in a low-g environment, (b) Time averaged heat transfer distribution in high-g (1.8 g) and low-g 
environments. 



Fig. 11. Schematic of axisymmetric liquid layer flowing on the inside of a tube (medium m). 
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All parameters in Eq. (19) can be measured with the exception of X/, 
from which the liquid film thickness 1/ can be determined. If heat 
transfer occurs across the liquid film, the temperature variation 
within the liquid needs to be accounted for. 


The range of film thicknesses that can be measured depends 
primarily on the absorption coefficient of the fluid. A plot of the 
effective transmissivity of a liquid film flowing within a silicon 
tube (Fig. 11) with 1 mm thick walls as the liquid film thickness 
is varied is shown in Fig. 12. The effective transmissivity is only 
about 0.33 even when no liquid flows in the tube due to reflection 
at the silicon interfaces. As the liquid film thickness increases, the 
effective transmissivity decreases and the system becomes effec- 
tively opaque at large thicknesses. If it is assumed that the trans- 
missivity can be accurately measured between 20% and 80% of 
these limits (horizontal dashed lines in Fig. 12), the range of liquid 
film thickness that can be measured can be found. For example, it 
is expected that only thin water films (~4-40 pm) can be mea- 
sured due to its high absorption coefficient while FC-72 films can 
be measured over a much broader range (~1 50-1 000 pm). 
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4. Closing remarks 

A new technique for the measurement of high resolution heat 
transfer distributions due to two phase flow phenomena based 
on IR thermometry has been described. The technique was demon- 
strated on droplet evaporation and flow boiling in low-g environ- 
ments. The equations needed to calculate the thickness of an 
isothermal liquid film in round tubes were derived and estimates 
were made of liquid film thicknesses that can be measured. 
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Appendix A 

A. t. Beer -Lambert law 

Consider radiation of intensity I 0 (W/m 2 ) striking a layer of 
thickness I. The intensity of the radiation (/) exiting the other side 
is given by J = / 0 exp(-aL) where a (m -1 ) is the absorption coeffi- 
cient and represents the ability of the layer to attenuate the incom- 
ing radiation. The transmissivity of a layer due to absorption alone 
can be expressed as t = exp(-aL). 

A.2. Emission from translucent materials 

Translucent materials simultaneously generate and absorb radi- 
ation. Radiation generated within the layer must travel through the 
material before it reaches the surface, during which it is attenuated 
in accordance with the Beer-Lambert law. Part of the radiation 
reaching the surface is emitted and the remainder can be reflected 
back into the layer where it can be re-absorbed and re-reflected. 
The effective radiant properties of translucent materials is derived 
below. It is assumed that the values all optical properties are inde- 
pendent of temperature. 

Consider first the radiation emitted by a differential layer of 
thickness dx with an emission coefficient s (m -1 ) and an absorption 
coefficient a (m _1 ) as shown in Fig. Al. The emission coefficient 
characterizes how black the layer is. The radiation flux emitted 
from the differential thickness dx in the x direction is given by 



Fig. 12. Range of liquid film thicknesses that can accurately be measured on a 
silicon tube. 


scr[T(x)] 4 dx. Because the differential element is thin, all of this 
radiation leaves the differential volume. The radiation that is not 
absorbed by the material that strikes the surface of the film at 
x = 0 is given by ea[T{x)] 4 exp[-otx]dx. The total radiation energy 
generated by the layer that strikes x = 0 can then be obtained by 
integration: 

E= I s(i[T{x)] 4 exp[-o(x]dx (A.1) 

Jo 

For the special case where the temperature in the layer is uni- 
form, Eq. (A.1) reduces to 

£ = ^(7T 4 [l -exp(-aL)] (A.2) 

In the limit of infinitely thick films, Eq. (A.2) indicates that the 
energy flux is 

E L ^=^aT 4 (A3) 

a 

Since very thick films act like a blackbody radiating energy at T, 
the layer will emit radiation according to El^oo = cl 4 , indicating 
that the relation between emission and absorption coefficients is 

s = ot. 

Consider next reflection at the surfaces. As mentioned above, 
part of the radiation reaching the surface of a layer can be 
reflected at the surface or transmitted through the surface multi- 
ple times as shown in Fig. A2 where medium m is sandwiched 
between two other media oo and n. The energy emitted by the 
layer that reaches the oo -m interface is given by E m and can be 
calculated from Eq. (A.1). The fraction of this energy that escapes 
from this surface is (1 p^^Em while the energy reflected back 
into m is p ix> _ m E m . The reflected energy that is not absorbed 
within m before arriving at the m-n interface is given by 
exp( -a m L m )p oo _ m E m = t mP^-rrfim- This radiation is reflected again 
at the m-n interface, absorbed within m, and arrives at the oo -m 
interface. The second contribution to the energy that escapes the 
oo -m interface is thus 

n - P»-m)Pco-mPm-„?m £ m (A.4) 

This process of reflection occurs again with the result that the 
third contribution is 

C* - Poo-m)Pi-mP 2 m-n'C E m (A.5) 

The energy escaping through the oo -m interface is therefore 


dx 

◄ — ►. 


a,e 


X 

►. 


L 

◄ ► 

Fig. Al. Radiation emitted within a translucent layer. 
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— P oo—m) P oo—mP m—n^ m Em ~~P * ‘ * 

OO 

= (1 — Poo-m^m 'y " , (P oo—mP m—n't m) 


(A. 6) 


Thus far, only the radiation emitted toward the oo-m interface 
has been considered. Radiation is also emitted toward the m-n 
interface and part of this energy escapes through the oo-m inter- 
face. By applying the same processes as shown in Fig. A2b, the en- 
ergy leaving the left surface can be expressed as. 

0 - P^-m)Pm-n T mEm + (1 - P oo-m) P m-n Poo-m T m E m + (1 

— Poo— m)Pm— nPoo— m^m^m + ‘ ' 


— 0 — P oo-m) P m-n 1 ' mEm ^ ( Poo-mPm-n^m ) 


(A. 7) 


The total emission from the oc-m interface (F m -oo) is the sum- 
mation of Eqs. (A.6) and (A.7): 

OO 

Em-o o - (1 - Poo-mK 1 + P m - n 1m)Em ^ {Poo-mPm-n^m}" (A.8) 

n= 0 

Since p m _ n , and z m are always less than unity, the infinite 
series may be written in analytical form as 


c _ 0 Poo— m )0 ’ Pm—n^m) n 

Zm-oo 1 _ 2 Zm 

r oo—mrm—n L m 


(A.9) 


Similarly, the energy emitted by m that escapes from the m-n 
interface (E m _ n ) is given by 


C _ 0 Pm— n)( 1 ~"P Poo—m^m) c 

c ra-n — 


(A. 10) 


1 - P oo—mP m—n^rn 

A3. Apparent reflectivity and apparent transmissivity of a single layer 

Let us first consider the reflection of incoming radiation from 
a single layer as shown in Fig. A3a. Part of the incoming energy is 
reflected from the oo-m interface. The remainder is attenuated by 
medium m, and part of this is reflected from the m-n interface back 
towards the oo-m interface. Using a similar line of reasoning as 
above to account the infinite number of reflections, the total 


energy reflected from m is the sum of the energies leaving the 
oo-m interface shown in Fig. A3a: 


Eoo-m — Poo- m Eoo + (1 - Poo— m) Pm-n T mE oo + (1 
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The apparent reflectivity is therefore given by 


(All) 


Poo—rr 


0 ~ Poo-m) Pm—nfJ 

1 Poo— mPm— n^m 


(A- 12) 


The apparent transmissivity of medium m to E ^ is found in a 
similar manner by summing the energies leaving the m-n interface 
in Fig. A3b: 


1 — Poo— mPm— n^m 


(A. 13) 


The apparent transmissivity is seen to be independent of the 
direction of the incoming radiation. 

The apparent absorptivity of m can also be determined from a 
similar accounting of energy absorbed within m to be 


0 - ftx,-m)d + Pm-n?m)(l ~ Tm) 

1 Pm—tlPo* — m^rn 


(A. 14) 


We leave it as an exercise to the reader to verify that 
0Capp,oo-m -EPapp,oo-m + Vp.oo-m = 1. These apparent properties are 
the same as those presented by McMahon (1950) if p oo _ m = p m _ n 
and the temperature distribution in medium m is uniform. 
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Fig. A2. Emission from medium m through the oo-m interface. Poo-m = reflectivity between medium oo and medium m, p m _ n = reflectivity between medium m and medium n, 
and T m = transmissivity of medium m. (a) Energy emitted by m towards the oo-m interface that leaves the oo-m interface, (b) Energy emitted by m towards the m-n interface 
that leaves the oo-m interface. 
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Fig. A3. Energy reflected and transmitted from medium m. (a) Apparent reflectivity, (b) Apparent transmissivity. 


A.4. Energy emission by two-layer materials including reflection from 
ambient 

Consider now the case shown in Fig. A4. The energy striking an 
IR camera is composed of four components: the energy emitted by 
an opaque surface s, emission from the two media (m and n), and 
by reflection from the surroundings. 

For the special case where s s = 1, we need to consider each of 
the four energy components below: 

AAA. Energy from surroundings 

The energy from the surroundings ( E ^ — crT^) is reflected from 
surfaces oo-m and m-n is given by p apP;00 _ m £oo- The energy trans- 
mitted through m is partly attenuated in n and the remainder is 
completely absorbed by surface s. No reflection occurs at surface 
s, simplifying the analysis. 


Ambient 
T ,£ =1.0 


Pm-n 


Surface 5 

£ s > T s=° 


Medium m 


Medium n 


Fig. A4. Schematic diagram for two-layer materials. 


A.4.2. Energy emitted by m 

The energy emitted by m leaving the interface oo-m is 
given by 


Em 



a m (7[T m (x)] 4 exp( 


a m x)dx 


(A. 15) 


The energy emitted by m that leaves interface m-n is again ab- 
sorbed either by n or s and does not reach the camera. 


A.4.3. Energy emitted by n 

The energy emitted by n that reaches interface m-n is 
given by 


E 


n 


a n a[T n (x)] 4 exp( -a n x)dx 


(A. 16) 


The energy that is transmitted through m and reaches the IR 
camera is given by z a pp,oo-mEn • 
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(A. 17) 


If the surface s is not black (e s < 1) but still opaque, reflections 
from this surface must be accounted for and the problem become 
more complex since an infinite number of reflections occur within 
two materials. For the sake of brevity, only the final results are gi- 
ven below. The energy striking the camera accounting for infinite 
reflections from surface s is given by 

Ec Poo-c^oo T~ Sm-cEm + Sn-cEn + ^s-c^s (A. 18) 


where 


AAA. Energy emitted by s 

The energy emitted from s ( E s — oT *) that reaches interface 
m-n is z n E s = exp(-ot n L n )E s . The portion transmitted through m 
that strikes the camera is z a pp,oo-m^nE s . 

The total energy striking the camera is the sum of these four 
components and is given by 
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Appendix B. Supplementary material 

Supplementary data associated with this article can be found, in 
the online version, at doi: 10 . 1016 /j.ijmultiphaseflow. 201 1 . 1 1 . 012 . 

References 

Coney, M.W.E., 1973. The theory and application of conductance probes for the 
measurement of liquid film thickness in two-phase flow. J. Phys. E: Sci. Instrum. 
6, 903. doi: 10.1 088/0022-3735/6/9/030. 

Dupont, V., Thome, J.R., Jacobi, A.M., 2004. Heat transfer model for evaporation in 
microchannels. Part II: comparison with the database. Int. J. Heat Mass Transfer 
47,3387-3401. 

Fischer, S., Herbert, S., Sielaff, A., Slomski, E.M., Stephan, P., Oechsner, M., 2011. 
Experimental Investigation of nucleate boiling on a thermal capacitive heater 
under variable gravity conditions. Microgravity Sci Technol. doi: 10.1 007/ 
si 221 7-01 1-9273-6. 

Gerardi, C., Buongiorno, J., Hu, L., McKrell, T., 2010. Study of bubble growth in water 
pool boiling through synchronized, infrared thermometry and high-speed 
video. Int. J. Heat Mass Transfer 53 (19-20), 4185-4192. 

Golobic, I., Petkovsek, J., Baselj, M., Papez, A., Kenning, D.B.R., 2009. Experimental 
determination of transient wall temperature distributions close to growing 
vapor bubbles. Heat Mass Transfer 45 (7), 857-866. doi:10.1007/s00231-007- 
0295-. 

Han, Y., Shikazono, N., 2009. Measurement of liquid film thickness in micro square 
channel. Int. J. Multiphase Flow 35 (10), 896-903. doi:10.1016/ 

j.ijmultiphaseflow.2009.06.006. 


Kattan, N., Thome, J.R., Favrat, D., 1998a. Flow boiling in horizontal tubes. Part 1: 
Development of a diabatic two-phase flow pattern map. J. Heat Transfer 120, 
140-147. 

Kattan, N., Thome, J.R., Favrat, D., 1998b. Flow boiling in horizontal tubes. Part 2: 
New heat transfer data for five refrigerants. J. Heat Transfer 120, 148-155. 

Kattan, N., Thome, J.R., Favrat, D., 1998c. Flow boiling in horizontal tubes. Part 3: 
Development of a new heat transfer model based on flow patterns. J. Heat 
Transfer 120, 156-165. 

Klausner, J.F., Zeng, L.Z., Bernhard, D.M., 1990. Development of a film thickness 
probe using capacitance for asymmetrical two-phase flow with heat addition. 
Review of Scientific Instruments 53(5), 3147-3152, 1992, doi: 10.1063/ 
1.1142568. 

Krebs, D., Narayanan, V., Liburdy, J., Pence, D., 2010. Spatially resolved wall 
temperature measurements during flow boiling in microchannels. Exp. Thermal 
Fluid Sci. 34, 434-445. doi:10.1016/j.expthermflusci.2009.05.005. 

Mani, P., Cardenas, R., Narayanan, V. 2012. Submerged Jet Impingement Boiling on a 
Nonuniformly Heated Polished Silicon Surface, Interpack201 1-52042, ASME 
2011 Pacific Rim Technical Conference and Exposition on Packaging and 
Integration of Electronic and Photonic Systems, Interpack 2011, Portland, OR, 
July 2011. 

McMahon, H.O., 1950. Thermal radiation from partially transparent reflecting 
bodies. J. Opt. Soc. Am. 40 (6), 376-380. 

Quinn, T.J., 1967. The calculation of the emissivity of cylindrical cavities giving near 
black-body radiation. Brit. J. Appl. Phys. 18, 1105-1113. doi:10.1088/0508- 
3443/18/8/310. 

Rouhani, Z., Axelsson, E., 1970. Calculation of void volume fraction in the subcooled 
and quality boiling regions. Int. J. Heat Mass Transfer 13, 383-393. 

Shedd, T.A., Newell, T.A., 1998. Automated optical liquid film thickness 
measurement method. Rev. Sci. Instrum. 69 (12), 4205-4213. 

Schweizer, N., Stephan, P., 2009. Experimental study of bubble behavior and local 
heat flux in pool boiling under variable gravity conditions. Multiphase Sci. 
Technol. 21 (4), 329-350. doi:10.1615/MultScienTechn.v21.i4.40. 

Sefiane, K., Moffat, J.R., Mater, O.K., Craster, R.V., 2008. Self-excited hydrothermal 
waves in evaporating sessile drops. Appl. Phys. Lett. 93, 074103. doi: 10.1 063/ 
1.2969072. 

Shen, J., Graber, C., Liburdy, J., Pence, D., Narayanan, V., 2010. Simultaneous droplet 
impingement dynamics and heat transfer on nano-structured surfaces. Exp. 
Thermal Fluid Sci. 34, 496-503. doi:10.1016/j.expthermflusci.2009.02.003. 

Thome, J.R., Dupont, V., Jacobi, A.M., 2004. Heat transfer model for evaporation in 
microchannels. Part I: presentation of the model. Int. J. Heat Mass Transfer 47, 
3375-3385. 


